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1. Introduction 



In spite of many recent developments, in particular in String Theory and Loop Quantum 
Gravity, the problem of quantization of General Relativity is still an open one both from 
mathematical and physical points of view. A key point of General Relativity is its close 
resemblance (and yet not equivalence) to an ordinary gauge theory. This is very clear 
in the first order (vierbein) formalism, where the dynamical variables are those suitable 
for a gauge theory of the Poincare group. This has been the starting point of the 
Loop Quantum Gravity approach too. This similarity led to various attempts to write 
down a regularized version of GR on a lattice. One of the most well-known is Regge 
Calculus PP, which is a natural and geometrically appealing way of discretizing gravity 
on a simplicial lattice. Recently it has been an active field of research in connection 
with Dynamical Triangulations, two dimentional Quantum Gravity, and therefore with 
conformal field theory and String Theory. One of the weak points of the simplicial 
lattice is that it is naturally related to the second order (metric) formalism of GR, so 
that its connection with gauge theory is rather obscure. The aim of this article is to 
continue the work started by M. Caselle, A. D'Adda and L. Magnea |2j, to reformulate 
Regge Calculus in terms of dynamical variables belonging to the Poincare group, so 
that it makes an explicit connection with gauge theory in general and Wilson-like gauge 
theory on a lattice jHj in particular. We have also found strict analogies between this 
approach and discrete Ashtekar variables jl]. We pursue this goal starting in section II 
by summarizing the main results of [2j in the second order formalism which consists of 
assigning the rotational degrees of freedom (Lorentz connection) to the links of a dual 
metric complex (Voronoi complex) of the original simplicial complex, so that they are 
functions of the translational degrees of freedom. In section III and IV we formulate a 
first-order principle, in the sense of Palatini f, in which the independent variables are 
the Lorentz connection and the normals to the n — 1-faces. These normals are considered 
analogues of the n-bein in the continuum theory, and the connection matrices as the 
connection one-form in GR. In section V we prove, in the case of small deficit angles, 
that Regge Calculus is a solution in the first-order formalism. This result is not obvious 
if we vary independently the two sets of the variables above. Moreover this equivalence 
for small deficit angles tells us that this first-order formalism has the effect of smoothing 
out some pathological configurations, like spikes, which affect Regge Calculus and might 
prevent the theory from having a smooth continuum limit. These spike configurations 
are in fact in the region of large deficit angles, where the first-order and second-order 
formalisms are not equivalent. In section VI we derive the general field equations for the 
connection matrices and for the n-bein. In section VII a measure for the path integral of 
this discrete theory of gravity is introduced and is shown to be locally invariant under 
SO(n). As a last step in section VII we propose a coupling of this discrete theory 

f The labeling of the first order method as Palatini has been questioned by F.W. Hehl, Four lectures 
on Poincare gauge field theory, Erice Cosmology Inst. 1979, pag. 40 footnote. See also Ferraris et al. 
...Palatini method discovered by Einstein 1925, Gen. Relat. and Grav. 14, pag. 243 (1982). 



Discrete Gravity in the first order formalism 



3 



of gravity with fermionic matter. This coupling is entirely performed by following the 
general prescription of the continuum theory. In other formulations of gravity (Regge 
Calculus and Dynamical Triangulations) the coupling with fermionic matter is usually 
introduced ad hoc. 

2. Second-order formalism 

In 1961 T. Regge [T] (see also [S] for a recent and updated summary of Regge Calculus 
and its alternative approaches) proposed a discretized version of General Relativity now 
known as Regge calculus. It is, mainly, based on the idea of substituting a Piecewise- 
Linear (PL)-manifold for differential manifold. In particular, these are simplicial 
manifolds built by gluing together two distinct n-dimensional flat simplexes by one 
and only one n — 1-dimensional face. The final product of this construction is called 
a simplicial complex, which owns the manifold structure if it has been made in such a 
way that each point of the simplicial complex has a neighbourhood homeomorphic to 
M. n JOj. The simplicial manifolds, we consider, are orientable. On each n — 2-dimensional 
simplex h, called a hinge, a deficit angle K(h) is defined 



where a l h , i — 1, p is one of the p n-dimensional simplices incident on h, and 9(o~h, h) 
its dihedral angles on h. The dihedral angle of a n-simplex on the hinge is the angle 
between the n — 1-dimentional faces that have the hinge in common. K(h) and the 
volume of the hinge V(h) can be expressed [Zj as functions of the squared length of 
the one-dimensional simplices (edges) of the complex. The Einstein-Regge action is, in 
analogy to the continuum functional over the simplicial manifolds and depends 

on the incidence matrix of the simplicial complex [Hj and on the (squared) lengths of 
the edges. Usually, in Regge calculus, the incidence matrix is fixed, so that the action 
can be written as 



The corresponding Einstein equations |Tj are derived by requiring that the action is 
stationary under the variation of the length of the edges. 

The original aim of this theory was to give approximate solutions of the Einstein 
equations in the case in which the topology is known. The theory is, as stressed by 
Regge himself, completely coordinate independent. 

It was pointed out in reference p j that a theory of Regge Calculus locally invariant under 
the Poincare group can be formulated by choosing in every simplex an orthonormal 
reference frame. In this way every n-dimensional simplex, considered as a piece of M. n , 
can be seen as local inertial reference frame, being flat, and an n-bein base can be chosen 
in it. As in the continuum theory, we need a connection between the reference frames 
of two simplexes that share a common n — 1-dimensional simplex. 




(1) 




(2) 



h 
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Figure 1. Dual Voronoi Plaquette 



Now we summarize and update some definitions and results of [2] that will be useful 
for our future discussions. Consider an hinge h and let {P^, ...,P^ l _ 1 } be its vertices. 
Suppose that this hinge is shared by k + 2 n-simplices (see FIG.l) {a, j3, 5\, 8%} whose 
vertices are labelled in this way 



a 



{Pi, P n -l, Qs^ai Qafd} 



(5 = {Pf, P n _ v Q a /3, Qj35\} 



h = { p i > -iPn-vQs^s^Qp^} ( 3 ) 

In each simplex S a we can choose an origin and an orthogonal reference frame. In this 
frame the vertices of the simplex S a have the following coordinates 



P = {Vi( a )} a=l,...,n , i = l,...,n-l 

Qs h k a = {^h a (a)} 

Qa/3 = {4f,(a)} (4) 

Given a n-dimensional simplicial complex, there exists a general procedure for defining 
the dual metric complex, called the Voronoi dual j2] |2I] [22] (see also (231 f° r an eas Y 
treatment and calculation at the details). The n- dimensional Voronoi polyhedron dual 
to a vertex P of the simplicial complex is the set of the points of the simplicial complex 
closer to P than to any other vertex, using the standard flat metric in the simplex. It 
turns out that the /c-simplex of the simplicial complex is dual to a (n — A;)-polyhedron 
in the dual Voronoi complex , and the A;-dimensional linear space identified by the k 
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simplex is orthogonal to the (n — k) dimensional space spanned by the corresponding 
polyhedron. In particular the point dual to the n-simplex S is the centre of a (d — 1)- 
dimensional sphere passing through all vertices of S (look at FIG.l). 
We can associate uniquely to a dual Voronoi edge an element of the Poincare group 
A(a,P) = {A£(a,/3), A a (a,f3)} by requiring that, 

A a b (a,(3)y b ([3)+A a (a,P) = y b (a) 

A a b (a,(3)z b ap (P)+A a (a,P) = z b a p(a) . (5) 

In other words we are embedding a and (3 in W 1 and adopting the standard notion of 
parallel displacement in M™. We move the origin of the reference frames of (3 to a. So 
that the position vectors in j3 of the vertices of the common face S a D Sp will become 
coincident with the position vectors of the same vertices in a. Thus the matrix A%(a, (3) 
will be an orthogonal matrix which describes the change from the reference frame of 
(3 to a, considered now as two different reference frames of the same vector space 
Since the simplicial manifold is orientable we can choose the reference frame in a and f3 
in such a way that A%(a,f3) is a matrix of SO(n). This procedure determines a unique 
connection that is the a torsion-free Levi Civita or Regge connection jH] • We can assign 
to each vertex D a , of the dual Voronoi complex, its coordinates x a (a) in the frame of 
a. Since we can choose arbitrarily the reference frames in a and j3, then the theory is 
invariant under arbitrary Poincare trasformations U(a), 

A(a,/3) i-> U(a)A(a, f3)U~ l (!3) 

x a {a) ^U£{a)x b {a) + U a {a) . (6) 

Now on we choose, following [2 J, to put the origin of each reference frame in the 
dual Voronoi vertices, then x a (a) = and the Poincare group is restricted to the 
rotation matrices without translations. Then A(a, (3) indicates only the rotation matrix. 
Without imposing this sort of gauge fixing, it might be possible to study a metric-affine 
theory of gravity on lattice and this could be the subject of a further research. See 
for more details on metric-affine theory of gravity well [Hj and [E] for the 

relative consequences in the continuum case. 

The hinge h is in one-to-one correspondence with its dual two-dimensional Voronoi 
plaquette that we still label by h. Now consider the following plaquette variable 

W a (h) = A(a,(3)A(P,5l)...A(5 h k ,a) . (7) 

As it has been shown [7j © is a rotation in the two-dimensional plane orthogonal 
to the n — 2-dimensional hyperplane spanned by the hinge h. The rotation angle is the 
deficit angle ((H). 

We consider n — 2 linear independent edge- vectors of the hinge h defined in the following 

way 

Et{a) = yt{a) - yU(a), E a n _ 2 {a) = y a n _ 2 {a) - y^a), (8) 
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and the antisymmetric tensor related to the oriented volume of h 

V Wab («) = ^^6^... Cn _ 2 Er(a)...K'r 2 2 (a) . (9) 

At this point it seems natural to propose the following gravitational action |2. 

/= -^E(V a2 {h) V [h) a ia2 (a)) • (10) 

h 

Let V(h) be the oriented volume of the hinge h. The action (jlUj) is equal [2] to 

I = ^sinK{h)V{h) , (11) 
h 

which for small deficit angles K(h) reduces to the Regge action (J2J). 
The presence in the action of sinK(h) instead of K(h) is only a lattice artifact 2 . The 
useful regime of this theory is sinK(h) pa K(h). In fact Frohlich jH] has conjectured, 
as claimed also in [2], that this action converges to the Einstein-Hilbert action, and 
J.Cheeger, W. Muller, R.Schrader in P have proved that the Einstein-Regge action (J2J) 
converges to the Einstein-Hilbert action. They proved that the convergence of the action 
is in sense of measure. This applies when, roughly speaking, the number of the hinges of 
the simplicial manifold increases along with the incident number of the simplices at each 
hinge. Then the triangulations will be finer and finer and the difference, in modulus, 
between the Einstein-Hilbert action on a manifold and the Regge-Einstein one, on the 
triangulations of the same manifold, become smaller and smaller. 

Let's consider the n — 1-dimensional face f a p = {Pi, ■■■,Q a /3} between the simplices a 
and (3. The vertex coordinates of f a p are xf(a), a = 1, n. The following vector [2j 



b a p „ (a) = W-fen-i^iW - x n (a)) bl ...{x n ^i{a) - x n (a)) bn 1 

= e o6l ... bn ^E bl ... E n n Jl (12) 

is normal to the face f a @, and it is assumed to point outward from the interior of the 
simplex a. Simple considerations of linear algebra show that b a p a (a) is an element of the 
inverse n x n, matrix built using the components of n linearly independent edge vectors 
of the simplex a multiplied by the determinant of the direct matrix (mathematically: the 
adjoint elements, modulo the determinant, of the matrix of the n-independent vectors). 
The analogous vector bp a a (f3) in the reference frame of j3 is related to the previous one 
by 

b af i a (a) = A a b (a,{3)b l3a b (f3) . (13) 
It can be easily proved that (JHJ) can be written as a bivector 

V(h) ClC2(a) = n\(n-\)\V(a) (^(a)^^^)-^^^^)), (14) 
where V(a) is the oriented volume of the simplex a. 
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We can consider the action as written on the dual Voronoi-complex of the original 
simplicial complex. The dual of the hinge h is a two dimensional plaquette whose 
vertices will be {a, /3, 8%, 8%j. These vertices are the dual of the n-simplices incident 
on the hinge h. 

If, briefly, we indicate A(ot, (3) as A a/ g, and so on, the holonomy matrix (JZj) around the 
plaquette is 

U aa = A af3 A f35 h...A s h a , (15) 

so that the action can be written in term of the bivectors b a p a (a), 

S = -± ][> (ff^VV)) ■ (16) 

h 

As it has been remarked in j2], this sum i) does not depend on the starting point of 
each two-dimensional plaquette chosen, and ii) it is a functional on the two dimensional 
plaquette of the Voronoi complex. 



3. First-order set up 

Up to now we have dealt with a second order formalism of discrete General Relativity 
. More precisely the connection matrices A a/ 3 and the b a p are both functions of the 
coordinates of the edges of the simplicial complex. Now we introduce a first order 
formalism in which A a p and b a p(a) are independent variables. We set, only, the following 
constraints 

b a Q0 (a)=A aP a b b b M > ( 17 ) 
which fix n independent conditions for each face, not enogh to determine the n ( n ~^ 
degrees of freedom of A a p. 

On the b a p(a) there is a further constraint since the n+1 normals to the n— 1-dimensional 
faces of a n-simplex a are linearly dependent, 

n+l 

^M") = ■ (18) 
0=1 

In this first-order formalism of General Relativity JT] we will consider torsion-free, 
metric-compatible connection matrices more general than the Levi-Civita one defined 
by equation (|5)i. (for a discussion on a possible weaking of the metricity condition see 
[TB] . Regge Calculus with Torsion has been studied by IT. Drummond [IT)] ) There is 
a technical problem we would like pointing out. The gravitational action in the form 
(fTT>j) could be dependent on the starting simplex . So we need to define the following 
antisymmetric tensor on the plaquette h 

W£l(a) = ^(V W («) + A a ^ h \(3)A Pa + ... 

+ A^...A^_ i ^vW(®A^ a2 _ i ...A fla ) c . (19) 




The action then is 

^^E Tr (W a )) ( 2 °) 

h 

The action (J2U)) coincides with the action ()16J1 in the second-order formalism, and, since 
(in matrix notation) 

WW(a) = A aP ..A Si _ lSi WW(5 i )A SiSi _ 1 ..Ap a . (21) 

equation (J2~U|) is independent of the starting simplex. 

Moreover the action (j2"Uj) is invariant under the following set of transformations 

^A^ = 0(a)A Q/3 0- 1 (/5) 
M«) ^ M a ) = 0{a)b aP {a) (22) 

where 0(a) and 0{(3) are two elements of SO{n). This can be interpreted as the 
diffeomorphism invariance of this discrete theory of Gravity. As we remarked before, 
without gauge-fixing, the affine group would have been the invariance group. This would 
have been the starting point to look for a description of Gravity, in the connection 
formalism, as a Yang-Mills theory of the affine group in the discrete case. We 
highlighted this problem as a, possible, subject for a future investigation. 

4. iV-bein on the simplicial complex 

It has been pointed out jHJEl that the h a p is the n-bein in the reference frame of each 
simplex. We are going to stress a feature which will be useful in the next. It has 
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been noticed in jH] that the way in which we have defined the coordinates, and, as 
consequence, the b a p itself, does not take into account the symmetry of each simplex. 
A better definition consists in introducing affine or "barycentric" coordinates. Each 
point P in the interior of a simplex (3 can be considered as a barycenter of n + 1 masses 
assigned at the vertices of the simplex (see reference ^0] and |Hj for more details). In 
particular, if we choose in the ro-simplex j3 a reference frame in which the origin coincides 
with the geometric barycenter [2|, the coordinates zf(/3) of the vertices, in this reference 
frame, by the definition of the geometric barycenter, satisfy the condition 

n 

J>«03) = O (23) 
which is analogous to (fT£|) . In barycentric coordinates [2j we may write the b a p as 

The meaning of this formula can be understood, geometrically, in an easy way if we look 
at FIG. 2. 63 ((]) is also normal to a face which belongs to 7. So it can be evaluated by 
solving (fTHj) . considered as an equation in 7, with respect to b 3 (/3). Equation (|2^j) does 
not give preference to the coordinates of any vertex. Thus the symmetry of the simplex 
is preserved as desired. The relation pij) can be inverted with respect to zf and, as in 
j2], we obtain 

= {n+ m nl V (a)r-i EC... aw - 1 4r i "- 1 ^(/^)---e- 1 1 (/3) • (25) 

It could be seen that (J23j) is related to the dual barycentric base as explained in [TU] 
and [Hj, but we do not explicitly discuss this link, since it is not relevant for the purpose 
of this article. 



5. First order field equations for small deficit angles 

As remarked in 0, in the second order formalism the action (|T6*j) is equivalent to the 
Regge action for small deficit angles K(h). In our first order formalism we don't have 
angles K(h). The variables related to the deficit angles are the connection matrices A a/ g. 
Then we assume, by definition, that the small deficit angle approximation in the first 
order formalism is the passage from the group variables h a p of SO(n) to the algebraic 
variables <p a p of so(n). As a consequence the connection matrices can be written in the 
form 

A a/3 = I + e <j) a p + o(e) (26) 

In order to avoid the technical complications, which we shall discuss in the next section, 
we now substitute the constraint b a p{a) = A a/3 bf3 a ((3) into the action for each Voronoi 
edge a(3, which, in the current approximation, is 

b aP (a) = (I + e <p a p)bp a ((3) + o(e) . (27) 
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A straightforward explanation of the action up to the first order shows that 

S = -i e J2 Tr + + <M °W h (a)) + o(e) (28) 

h 

where °W C1C2 is the (2, 0) antisymmetric tensor at zeroth order, in the e- 

expansion. Here we have used the approximation ()26|) . For each Voronoi-edge we have 
solved the constraint (fTTj) up to the first order (|2*7jl . 

We also impose the requirement that the action is stationary under variation with respect 

tO <p af3 , 

8S 



e 



J2 °W h C1C2 + o(e) = . (29) 



These are the implicit equations for (p a/3 , which has nt ^ n ~ 1 ^ independent 

components. Since 



d°W h C1, » (<!>„*, bi 



b aP I n\(n-2)\V(P) 



( S Cl h Cl h s — A Cl h C2 h 



— l^h Cl h s — ft C2 h Cl h r 
i 0„ u /3Sh U aP °s u p5h u af3 



(30) 



the determinant of this matrix is non zero (notice that matrix (jSUj) has dimension " ^ ' ). 
By inverse function theorem (also labeled Dini's theorem), locally, we can invert equation 
f!29|) . so that, again locally, there is one solution giving (J) a p as a function of the n-bein 
6j. We do not know anything about the uniqueness of this solution globally. We are, 
in some sense, in a situation analogous to that described in [T2]. Barrett found, locally, 
the uniqueness of the solution, but has showed that this result is not valid globally. 
Locally the unique solution can be found by using the Levi Civita — Regge connection. 
We determine the matrix connection by equations (0), in which the translation vector 
is a known quantity, once we put the origins of the refernce frames in the barycenter of 
the simplices. The equations (0) determine the connection matrices as functions of the 
coordinates of the vertices. In barycentric coordinates these are in one-to-one relation 
with the hi. Then we determine the connection matrices as functions of the hi. The 
identity 

V h ClC2 {a) = , (31) 

he(aP) 

can be easily proved if we express 

V (h) C 1C2 

(a) as in (fllj) . and use the constraint 
Y^B=x = 0- The statement that this is the Levi Civita- Regge connection implies 

V< & >(«) = A a ^ h \(3)Apa = = A^...A 5ti5 ,V^(^)A^ ti ...A^ , (32) 
so that to the zero order 

V ih \a) = V {h) {(3) + 0(e) = ... = V {h \5 h k ) + O(e) . (33) 
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These facts imply that 

e ^ °W h ClC2 = ( eVk ClC2 (°) + e °( e )) = + °( e2 ) ' ( 34 ) 

that is to say the Levi Civita-Regge connection is the solution of our first-order equations 
for the connection matrices in the limit of small deficit angles. 

6. First order field equation: the general case 

In the previous section we have seen that, in the case of small deficit angles, Regge 
Calculus is the solution of the first order field equations. 

Now we are going to deal with the general problem. We would like to derive the equation 
of motion by varying the action with respect to A a p and b a p. First we have to take in 
account the constraints (fTTjl and (JTBJ). So, in order to perform independent variations of 
A a /3 and b a p, it is necessary to put Lagrange multipliers in the action. Then the action, 
in the Palatini first order, is 



h (a/3) 

+ ]TTr {\^(A af sA T a(S - /)) + $>(«) (EM«) ] , (35) 

(a/3) a \/3=l / 

where A^- 1 and fj-(a) are n-dimensional vectors and Lagrange multipliers respectively, 
and \( a/3 *) is an n x n matrix. The constraint A aj gA^ — I is introduced to restrict the 
variation of k a p on the group SO{n). 
We introduce the following lemma: 
Lemma: The action, in matrix notation, 

S' = Tr(AA) + Tr (A(AA T - /)) (36) 

gives the equation of motion, (if we assume that the variation with respect to A is 
stationary) 

(AA) = (AA) T . (37) 

Proof: If we consider the variation of the action, using the property that Tr(M) = 
Tr(M T ) and the action is stationary respect to A and A T we have 



= ^ = A + A T X , 
5A 

= TrF = A T + XA . (38) 
oA 1 

Multiplying the first equation for A on the left, and the second for A T on the right, and 
subtracting term by term these two equations, we have (}3Tj) . 
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Applying the lemma to the action (|35|) for the variation with respect to A a ^, we obtain 
the field equations 

E (ULW^a))^ - \ aP t b aP {a) j = J2 (U h aa W h {a))l - X af3 ,M«) , . (39) 

he (a/3) he(a/3) 



The next step will be to determine the field equations for the variations of b a p. For 

. As remarked in ;2|, we h 



this it is necessary to determine the quantity 9 g b a ■ As remarked in [2|, we have the 
following identity 



(Via))- 1 = ^6 a ,.. 0n e^-^6-(/3)...6-(^) . (40) 

We can write this formula in a way which is not dependent on the chosen index 
j.Equivalently equation ()4()j) can be written as 

1 1 n+1 



(Via))"" 1 = —-—^e^ e^%{f3)...b%{(3) , (41) 
n! n + 1 ^— ' J J 

so that it is clearly independent of any index. Thus, we have 



Via) 



1 1 n+1 



Evaluating its derivative with respect to bf(a) we obtain 



86?(a) n\in - l)V 2 ia) aK ' ' 1 J 

Now we are ready to derive the field equations for the variations of the b a p. So let's 
consider the action (|35|). It is straightforward to write it as 



S = ~\j2 Tr K«W h ia) + U T h aa W T \a)) + E (M«) - A^ft/toGS)) 

ft (a/3) 

+E^ (a M) (a q/3 a^ - /)) +£>(«) Em«) • ( 44 ) 

(a/3) a \/3=l / 

Let's define 

the variation with respect to b l aj3 ia) is 

= E 4n\in-2)\V(af< ij ~ n\(n - l)V {af™ ^^)^Yj 
+\ a/3 - m(a) = . (46) 
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The field equations (pftj|) and (|3U|) are the analogues of the first order Einstein equations 
in the n-bein formalism (see [IS] for a strict analogy). The connection matrices, A aj g, 
the n-bein, b l a Ja), and the Lagrange multipliers X l a/3 and (J? (a) must be determined by 
solving the equations (j4T)j) and (J3l?j) and the constraints (fTTj) . (fTHj) as well the orthogonal 
condition of the connection matrices. 

A nice feature of these equations is that if |A a/ g, b l a {(3), \ a/ 3,/ij(a;)| is solution of the 
equations of motion, the map 



\ a/3 l— * ^ Aj a/3 

Wi(a) h-> uT^/i^a) . (47) 

will provide an other solution, as it is easily to verify, of the discrete Einstein equations. 
This map, if we restrict to the b's only, shows that the discrete Einstein equations have 
a conformal symmetry, anyway in the solutions not in the action. This is the same 
feature of the Regge equations pQ. 

7. The measure 

In this section we shall discuss the quantum measure to associate with the previous 
classical action. The action [2] is invariant under the action of the group SO(n), that 



is the gauge group 
We will use the following notation: 

Ai(M«)) = dbl p {a)...dbl p (a) , (48) 

and let 

w(A Q/3 ) (49) 
be the Haar measure on SO(n). The partition function for this theory is: 

/n+1 
e l Eh Tr {u^ a w h {a)) ^(^^ b a p(a)) J^J 5{b a p(a) — A a/ g6 i a a ( ( g))«(A a/ g)/i(6 Q , / 3) .(50) 

a 13=1 a/3 

Here the product JT is a product over all vertices of the dual complex as well the product 
W a a over all the Voronoi links which has a as one of his vertices. It is straightforward 
to see that the measure is invariant under the gauge transformation ()22j). In fact if we 
perform the gauge transformation (|22j), the modifications to the measure are: 



n+1 

5(X; & a/3(«)) 5 (^(«)-AU^(«)^(Aa/3))"(^) ■ (51) 
(3=\ 
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But n{b' a p) is equal to det\0(a)j n(b a p), and reduces to ^{b a p) as well. The Haar 
measure of SO(n) is right and left invariant. So fi(A' a/3 ) = yu(A Q/3 ). Again, by using 
equations (j22j) the two deltas of (f5*T]) can be written in the following form: 



tf(0(a)£>«,(a))) , 
p=i 

6[0(a)(b al3 (a) - Aapbpa^))) . (52) 

Then the properties of the delta function along with the last considerations prove that 
equation (JoTj) can be written as 

n+l 

S(^b a p(a))5(b a p(a) - A a/3 bp a ^))fj,(A a p)ii(b a p) , (53) 

13=1 

This establishes the invariance of the measure under gauge trasformations. 
Challenged by [TH] we can argue that the metric structure we are considering is not as 
peculiar as the metric written as a function of the edge-lengths in Regge Calculus. In 
first-order formalism we neither have transition functions which depend on the metric 
structure, the n-bein in our case, nor do we sum over a metric that is gauge-fixed. From 
the form of the constraint equations, those in the argument of the delta functions, and 
from the calculations we have performed, we expect that this measure, once we have 
integrated over the deltas, is really highly not local as has been discovered for Regge 
Calculus. 



8. Coupling with matter 

In the continuum theory on Riemannian manifolds with torsion-free connection the 
coupling with fermionic matter (for the same case in presence of torsion see ^3]) is 
given by the lagrangian density 

C = % - (V<7 a V^ - <V^7» - rnW • (54) 

The 7 a , a = 1, ...n are the Dirac- matrices satisfying the Clifford algebra 

7 <y + 7 y = 25 ab , (55) 

whereas if) is the n-dimensional Dirac spinor field = ^^7 1 ),V At the covariant derivative, 
and the n-beins on the tangent space of the Riemannian manifold (M,g), where the 
lagrangian density is defined (jM}, such that 

gr(x) = eZ{x)4{x)8* . (56) 

We are assuming that the Riemannian manifold (M, g) in question has a spin structure, 
that is its second Stiefel- Whitney class is zero. 
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Now we have all the ingredients to define the coupling of gravity with fermionic matter 
on the lattice in analogy with the continuum case. Let 2v = n or n = 2u + 1 (depending 
on whether n is even or odd) and consider the 2 1/ -dimensional representation of the 
two-fold covering group of 5'0(n)jH]. So, instead of considering the connection matrices 
A Wj g, we will deal with the 2 U x 2 U connection matrices D a p such that 

DarfD^ = (Aa ) a b j b . (57) 

Given D a p we can determine A a/ g. Furthermore if we know A Q/ g, we can determine D a p 
up to a sign. In particular, from (|57)l . we can write A a p as [E] 

{Kp) ah = ^Tr( la D aPlb D-D (58) 

In the discrete theory we assume that the spinor field is a 2 V complex vector defined at 
each vertex of the dual Voronoi complex, that is to say a map that to each vertex a 
associates the 2 V complex vector ip(a). 

In order to define the covariant derivative on a lattice we have to derive the distance 
| a/3 1 between the two neighboring circumcenters in a and f3. Our reasoning concerning 
the baricenters can be extended to the circumcenters too. The distance Ah\ of the 
circumcenter in a from the face afl can be determined by calculating the volume of the 
n-dimensional simplex obtained by joining the circumcenter with the n-vertices of the 
face a/3, and dividing it by n and by the volume of the face itself. So we have 

n A \o a0 {a)\ 

in which the are the circumcentric coordinates of the vertices of the face a/3 and 
|& a/ g(a)| the module of the b a p(a) written as function of the zf, as in 
In the same manner we have 

Afl2 = ~3 \h (aw ' 60 

Finally, we have \a/3\ = Ah\ + A/12. 

We are ready to define the covariant derivative (V A1 '0)(a) on a lattice 

(V^)(a) = ^ . (61) 

So far the discrete version of the the action for the coupling between gravity and 
fermionic matter can be written in the form (see also [2*U] ) 



v (a/3),/3=l,...,n+l 

-D a/3 ^(P)b a a ^Ma)^) - mV(a)tP(a)iP(a)^ . (62) 
Then the quantum measure, which also includes fermionic matter, can be written as 
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/n+l 

a /?=1 

x Yl 5(b af3 (a) - A a pbp a ^))n(D af3 )n(b a p) (63) 

a/3 

where S is the action for pure gravity (J20J), (j,(D a p) the Haar measure on the two-fold 
covering group of SO(n), while n(ip(a)) = dif)(a) is the standard measure on C T ' . 

9. Conclusions 

In this paper we have studied a discrete theory of gravity in its first-order formalism. 
Following an earlier previous paper (2] , we have chosen an orthonormal reference frame in 
each simplex and have defined a connection matrix as a transformation matrix between 
the reference frames of two n-dimensional simplices that share a common n — 1-face, 
considered as two distinct reference frames at the same point. This defines a Levi Civita- 
Regge connection. These matrices allow us to write an holonomy matrix around each 
hinge h. We define the action as the sum over the hinges of the traces of the holonomy 
matrices multiplied by the oriented volumes of the hinges. This action can be written 
on the dual Voronoi complex of the original simplicial complex. We can express the 
action as a function of the connection matrices and of the vectors b a p which are normal 
to the faces of the simplices and whose modulo is proportional to the volume of the 
faces themself. 

The action is very similar to the Wilson action of lattice gauge theory. Here the b a p 
have the same role as the n-bein in the continuum theory. Moreover the action is locally 
invariant under the action of the Poincare group. 

On the dual Voronoi metric complex the first-order formalism is implemented by 
considering the dynamical variables A a/ g and b al 3 independent. It is shown that in 
the limit of small deficit angles the Levi-Civita or Regge connection is a solution of the 
equations of motion. As in Barrett [12] . the solution is unique locally, by Dini's theorem. 
However we do not know what is going to happen globally. 

The general equations for A a/3 and b a p are derived by using Lagrange multipliers, in 
order to take in account the constraints of the first order formalism. 
A quantum measure and the relative partition function has been defined. They are 
locally invariant under the action of SO(n). We have introduced a coupling of gravity 
with fermionic matter on the Voronoi complex as well. This coupling seems as natural 
as for the continuum theory. In particular it seems the Voronoi complex allows us to 
introduce the coupling with matter, avoiding all the troubles we had in Regge Calculus. 
We hope that this first-order formalism might be useful for numerical simulations. In 
particular its first-order character, numerically, could have more advantages than the 
usual (Regge) second order formulation in implementing the evolution of equations of 
motion. 
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